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1 Introduction 

It is well known that the classical Black-Scholes model does not allow in real applications to replicate perfectly 
contingent claims. Of course, this is due to market incompleteness and specifically two major reasons : the 
non-Gaussianity of prices log-returns and the finite number of trading dates. The impact of these features 
have been intensively studied separately in the literature. 

There is a large literature on pricing and hedging with non Gaussian models (allowing for stochastic 
volatility or jumps), in a continuous time setup. Then, the hedging error related to the discretization of 
the hedging strategy is in general ignored or investigated separately. One popular approach is the Variance- 
Optimal hedging. Let S c denotes the underlying price process where the superscript c refers to the continuous 
time setting) ; if H denotes the payoff of the option, the goal is to minimize the mean squared hedging error 



over all initial endowments c£l and all (in some sense) admissible strategies v. The first paper specifically on 
this subject is due to Duffle and Richardson, see [T8]. Among significant early contributions there are [31)1 [571 
I39U331[2"5] . a fairly complete recent article on the structure of mean- variance hedging, with a rich bibliography 
is provided by [TTj. One of the now classical tools is the so called Follmer-Schweizer decomposition. Given 
a square integrable r.v. H and an (J-t)-semimartingale S = (St)t>o, that decomposition consists in finding 
a triple (Hq^,L) where Hq is Jxi-measurable, £ is (J r t)-predictable and L is a martingale being orthogonal 
to the martingale part M of S such that H = Hq + J Q T £ s dS s + Lj>. In the recent years, some attention was 
focused on finding explicit or quasi explicit formulae for the Follmer-Schweizer decomposition or the optimal 
strategy for the mean- variance hedging problem. For instance [5] gave an expression based on Clark-Ocone 
type decompositions related to Levy type measures when the underlying is a Levy martingale, |15| still in 
the martingale case with techniques of partial integro differential equations. [25] obtained significant explicit 
decompositions when the underlying is the exponential of a Levy process and the contingent claim is a vanilla 
type option appearing as some generalized Laplace transform of a finite complex measure. Other significant 
semi-explicit formulae appear in |301 [31] . [29] was continued by [28] in the framework of processes with 
independent increments with some applications to the electricity market. 

However, in practice, the hedging strategy cannot be implemented continuously and the resulting optimal 
strategy has to be discretized. Hence, to be really relevant the hedging error should take into account this 
further approximation. 

An alternative approach, less investigated in the literature, is to consider directly the hedging problem in 
discrete time as proposed by Cox Ross and Rubinstein [16] . The first incomplete market analysis in the 
spirit of minimizing a quadratic risk is due to |19| . They worked with the so-called local risk- minimization. 
The problem of Variance-Optimal hedging in the discrete time setup was proposed in |35l I38| . In the recent 
years some interest on discrete time was rediscovered in [5J [HI [32] • [12] revisits the seminal paper [TH] in the 
spirit of global risk minimization. In the discrete-time context, a significant role was played by the analogous 
of the previously mentioned FS-decomposition. It is recalled in Definition 12.81 

Recently, many approaches have been proposed to obtain explicit or quasi-explicit formulae for computing 
both the variance optimal trading strategies and hedging errors in discrete time. For instance, in [T], Angelini 
and Herzel derive closed formulae for the variance optimal hedge ratio and the corresponding hedging error 
variance when the underlying asset is a geometric Brownian motion which is martingale. As we said, Kallsen 
and co-authors contributed at providing semi-explicit formulae for the Variance-Optimal hedging problem 
both in discrete and continuous time, for various kind of models. In particular in |29| . semi-explicit formula 
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are derived for the (discrete and continuous time) Variance-Optimal hedging strategy and for the resulting 
hedging error, in the specific case where the logarithm of the underlying price is a process with stationary 
independent increments. One major idea proposed in [29] and |10| consists in expressing the payoff as a 
linear combination of exponential payoffs for which the variance optimal hedging strategy can be expressed 
explicitly. With a similar methodology and in the same setting, Angelini and Herzel [3] determine the 
Laplace transform of the variance of the error produced by a standard delta hedging strategy when applied 
to several class of models. In [IT] similar results are provided in the continuous time setup. In this paper, 
we use the generalized Laplace transform approach to extend the results of [53] to the case of processes with 
independent increments (PII) relaxing the stationary assumption on log-returns. The semi-explicit discrete 
Follmer-Schweizer decomposition is stated in Proposition 13. Ill the solution to the mean- variance hedging 
problem in Theorem 14.11 The expression of the quadratic hedging error in Theorem 14.31 gives a priori a 
criterion of market completeness as far as vanilla options are concerned. This confirms that the (even not 
stationary) binomial model is complete, see Proposition 14.51 

Our discrete time model consists in fact in the discretization of continuous time models which are exponen- 
tials of processes of independent increments. Given a continuous-time model (S' t c )t>o, where S% — sq cxp(V t c ) 
and X c is a process with independent increments and discrete trading dates to, ti, • ■ • , tjsr, our discrete model 
will be S = (Sk), such that Sk = S% , for all k = 0, 1, • • • N. In this discrete time setting, the Variance- 
Optimal pricing and hedging problem consists in looking for the initial endowments c £ R and the admissible 
strategy v = («&) which minimizes 

N 

E[(Vfi - Hf] with V$ = c + v k AS k . 

fe=i 

This framework is indeed well suited to take into account together both the non-Gaussianity of log-returns 
and hedging errors due to the discreteness of trading times. Our investigation for quasi-explicit formulae 
when the underlying is the exponential of sums of independent random variables is due to two reasons. 

1. The first one comes from the fact that the basic continuous time model can be time-inhomogeneous in 
a natural way, see for instance |28) . 

2. The second, more original reason, is that the discretized times, which correspond in our case to the 
rebalancing dates, are not necessarily uniformly chosen. 

About item 1., some prices exhibit non stationary and non-Gaussian log-returns. One common example of 
this phenomenon can be observed on electricity futures or forward market: the forward volatility increases 
when the time to delivery decreases whereas the tails of log-returns distribution get heavier resulting in huge 
spikes on the Spot. The exponential Levy factor model, proposed in [7] and [13] allows to represent both the 
volatility term structure and the spikes on the short term. More precisely, the forward price given at time t 
for delivery of lMWh at time Td > t, denoted Fj d is then modeled by a two factors model, such that 

S r t := F?* = F^ d exp(mp + / a s e~ x ^~ u UA u + a L W t ) , for all t e [0, T d ] , (1.1) 

Jo 

where m is a real deterministic trend, A a real Levy process and W a real Brownian motion. Hence, forward 
prices are modeled as exponentials of PII with non- stationary increments and existing results from [29J valid 
for stationary independent processes cannot be applied for that kind of models. 

Concerning item 2., the announced motivation for our development is to be able to analyze the impact of 
a non-homogeneous discretization of the trading dates on the Variance-Optimal hedging error. The issue of 
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considering non-homogeneous trading dates was first considered by Geiss [21] who analyzed the impact on 
the hedging error of discretizing a continuously rebalanced hedging portfolio. He showed that for a given 
irregular payoff (e.g. a digital call), concentrating rebalancing dates near the maturity instead of rebalancing 
regularly can improve the convergence rate of the hedging error. Later, Geiss and Geiss [55] introduced the 
so called fractional smoothness quantifying the impact of the payoff irregularity on the optimal discretization 
grid. The reader can consult [23] for a nice survey on this subject and [23] for some recent developments. 

Hence, it seems to be of real interest to be able to consider such non-homogeneous grids. However, if 
the continuous time log-price model X c = log(5 c ) — log(so) na s independent and stationary increments, 
considering non-homogeneous trading dates involves a non stationary discrete time process X such that 
Xk = X ( { k for k = Q,---N, where to,ti,-- - , tjy denote the non- homogeneous trading dates. Hence, here 
again existing results from |29) cannot be applied neither for hedging at non-homogeneous times nor for 
evaluating the resulting hedging error. 

In the present work, we have performed some numerical tests concerning both applications. One major 
observation is the remarkable robustness of the Black-Scholes strategy that still achieves quasi-minimal 
hedging errors variances, with both non Gaussian log-returns and discrete rebalancing dates. Besides, our 
tests show that when hedging with electricity forward contracts, the impact of the choice of the rebalancing 
dates on the hedging error seems to be more important than the choice of log-returns distribution (Gaussian 
or Normal Inverse Gaussian, in our case) . Concerning the case of hedging an irregular payoff (a digital call, 
in our case), our numerical tests confirm the result of |21| . In almost Gaussian cases, we observe that the 
variance optimal hedging error, can be noticeably reduced by optimizing the rebalancing dates. However, 
this phenomena is less pronounced when the tails of the log-returns distribution get heavier for which the 
hedging error gets less sensitive to the rebalancing grid. This suggests that the result of [5T] and [M] could 
not be extended straightforwardly to the non Gaussian case. 

This article is organized as follows. In Section 2, notations and generalities on the discrete Follmer- 
Schweizer decomposition are presented. In Section 3, we derive semi-explicit Folhner-Schweizer decom- 
position for exponential of PII. Section 4 is devoted to the solution to the global minimization problem. 
Illustrative example and simulation results are given in Section 5; in particular, subsection 15.21 is concerned 
with data coming from the electricity market. 

2 Generalities and Discrete Follmer-Schweizer decomposition 

We present the context of the problem studied by [55] . 

Let (£1,7-, P) be a probability space, N £ N* a fixed natural number and F = (J-k)k=o,--- .n a fixed 
reference filtration. We shall assume that T — Fn- Let (Sk)k=o,--- ,n be a real-valued, F-adapted, square- 
integrable process. We denote by ASk the increments Sk — Sk-i, for fc = 1, • ■ ■ , N. We use the convention 
that a sum (respectively product) over an empty set is zero (resp. one). 

Definition 2.1. We denote by O the set of all predictable processes v (i.e.: vj~ is J- k-i -measurable for each 
fc > 1) such that WfcAS'fc 6 £ 2 (Q) for fc = 1, • • ■ , N. For v £ Q, G(v) is the process defined by 

k 

G k (v) := ^UjASj, for k = l,---,N. 

i=i 

The problem addressed in [38] is the following. 
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Given H G £ 2 (f2), we look for (Vg*, (p*) which minimize the quantity 

e[(H -V -Gt&)) 2 ] , (2.2) 

over Vg G R and (p G 0. It will be called discrete time optimization problem. The expression 
E 



will be called the variance optimal hedging error. 



Definition 2.2. Schweizer ]°J8^ introduces the following non- degeneracy condition (ND). We say that 
S satisfies the non- degeneracy condition (ND) if there exists a constant 8 G]0, 1[ such that 

[E[AS k \F k ^]) 2 < 8E[{AS k ) 2 \T k ^] , 

P.a.s for k = 1, • • • ,N. 

Remark 2.3. 1. If (S k ) is a martingale then (ND) is always verified. 

2. Note that by Jensen's inequality, we always have (E[A5fc |J-fc_i]) 2 < E[(AS , fe) 2 |J-/ £ _i] a.s. The point of 
condition (ND) is to ensure a strict inequality uniformly in uj. 

To obtain another formulation of (ND), we now express S in its Doob decomposition as S k = M k + A k 
where M k is a square-integrable martingale and A k is a square-integrable predictable process with A® = 0. 
It is well-known that this decomposition is unique and is given through 

A^Lfe := E[A5 fe |J^_i] , and AM U := AS k - AA k . 
We will operate with the help of some conditional moments and conditional variance setting 

Var[AS k \T k -i] := E[(AS k ) 2 \ F k ^] - E[AS k \ T k -i] 2 . 
Remark 2.4. For k = 1, . . . , N , we have the following. 

1. E[(AS k ) 2 \F k ^] = E[(AM k ) 2 \T k ^} + (AA k ) 2 ; 

2. Var[AS k \F k -x] = E[{AM k f\F k _ 1 ] ; 

3. Previous conditional variance vanishes if and only if AM k = a.s. 
We introduce the predictable process X k by 

AA k E[AS k \T k ^] 
k ■ E[(A^) 2 |^-i] E[(AS , fe) 2 |J r fe _i] ' 1 ^ 

for all k = 1, ■ • • , N. These quantities could be theoretically infinite. 

Remark 2.5. Suppose that P(AS k = 0) = for any k = 1, • • • , N. 

1. Then E^ASk) 2 ^^] > a.s. In fact, let B = {u\E[(AS k ) 2 (w)\ F k -i] = 0}. This implies AA k = 
on B because of Remark \2.4\ 1- By the same Remark, 

= l B E[(AM k ) 2 \T k ^] = E[l B (AM k ) 2 \T k ^] , 

so AM k = a.s. on B. This implies that AS k = a.s. on B. By assumption, B is forced to be a null 
set. 
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2. Previous point 1. guarantees in particular that (Afe) are all finite. 



Definition 2.6. The mean-variance tradeoff process of S is defined by 



j 



E[AS l \T l . 1 }' 2 

Var[ASi\Fi-i] ' 



for all j = 1, • • • ,N. K d is the discrete version of the continuous time corresponding process K defined for 
instance in Definition 2.11 of ]28^ or in Section 1. of ]3Gj . 

Proposition 2.7. The condition (ND) is fulfilled if and only if 



A basic tool for solving the optimization problem (|2.2p in [38] is the discrete Follmer-Schweizer decom- 
position. 

Definition 2.8. Denote by S = M + A the Doob decomposition of S into a martingale M and a predictable 
process A. A complex-valued square integrable random variable H is said to admit a discrete Follmer- 
Schweizer decomposition (or simply discrete FS- decomposition) if there exists a To-measurable Hq 7 a 
complex-valued process £ such that both Re£(z), Im^(z) belong to 0, and a square integrable C-valued mar- 
tingale L H such that 

1. L H M is a martingale; 



When Point 1. is fulfilled L and M are called strongly orthogonal. 

If H is a real valued r.v. then H admits a real discrete FS decomposition if it admits a FS decomposition 
with Hq G R and £ being a real valued process. In this case £ £ O. 

2.1 Existence and structure of an optimal strategy 
Assumption 1. (Sk)k=i,--- ,N satisfies the non- degeneracy condition (ND). 

Remark 2.9. 1. Under Assumption [7J Proposition 2.6 of 1381 guarantees that every square integrable 
real random variable H admits a real discrete FS- decomposition. 

2. That decomposition is unique because of Remark 4- 11 of \33j . 

3. The previous two points imply the existence and uniqueness of the discrete Follmer-Schweizer decom- 
position when H is a complex square integrable random variable. 

4- An immediate consequence is that the decomposition of a real square integrable random variable is 
necessarily real. 



Vor[A5 fc |Jfc-i] 



is a.s. bounded uniformly in uj and k. 



Proof. See (1.6) in [Ml- 



□ 



2. E{L$) = 7 

3. H = H Q + £ fe A5 fe + L 
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Other tools for solving the optimization problem and evaluating the error are the following. 
Proposition 2.10. If S satisfies (ND), then Gjv(0) is closed in £ 2 (P). 

Proof. See [3S], Theorem 2.1. □ 

Theorem 2.11. Suppose that S = M + A has a deterministic mean-variance tradeoff process. Let H be a 
square integrable real random variable with discrete real FS- decomposition given by H = Hq + Gn(^ H ) + L^. 

1. The optimization problem 112. 2\) is solved by (V *,ip*) where V * = H and ip* is determined by 

<Pk = $ + A fc (F fc _i -Ho- G fc -i (*>*)). 

2. Suppose that is a trivial a -field. The hedging error is given by 

N N 

fe=l j=k+l 

Proof. Point 1. follows from Proposition 4.3 of [35]. Concerning Point 2., Lq = a.s. since Tq is trivial. 

The result follows from Theorem 4.4 of [35]; □ 

Similarly to [22], we will calculate it explicitely in the case where S is the exponential of process with 
independent increments. 

3 Exponential of PII processes 

From now on, we will suppose that (X„)„ = o i ... ; jv is a sequence of random variables with independent 
increments, i.e. (Xi — Xq, ■ ■ • ,Xn — Xjv-i) are independent random variables. From now on, without 
restriction of generality, it will not be restrictive to suppose Xq = 0. We also define the process (S n ) n =o,--- n 
as S n = sq cxp(X n ), < n < N for some so > 0. 

Definition 3.1. We denote D = {z e C\ exp(zX N ) e C 1 }. 

3.1 Discrete cumulant generating function 

Definition 3.2. We define the discrete cumulant generating function as 

m : D x {0, • • • , N} — > C with m(z, n) = K[e zAXrl ] for all n = 1, • • • , N and by convention m(z, 0) = 1. 

This function is a discrete version of the cumulant generating function investigated in |28) . 

Remark 3.3. 1. If z g D then the property of independent increments implies thatm(z,ri) = E[cxp(zAX„)] 
is well-defined for all z G D and n = 0, 1, ■ • • , N. 

2. If j 6 R + n£>, Cauchy-Schwarz inequality implies that [0, 7]+* C D; if 7 e R~CiD then [7, 0]+iM. C D. 
This shows in particular that D is convex. 

Remark 3.4. When X has stationary increments then we have m(z, n) = m(z, 1) for all n = 1, • • • , N . We 
denote this quantity by m(z) similarly as in \29$ . Section 2. 

We formulate some assumptions which are analogous to those in continuous time case, see [2"5] . 
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Assumption 2. 1. AX n is never deterministic for every n = 1, • • • , N. 
2. 2eD. 

Remark 3.5. In particular, S n G C 2 (fl), for every n = 0, 1, • • • , N , because 2 G D. 

Lemma 3.6. z i— > m(z, n) is continuous for any n = 0, 1, • • ■ , iV. particular, if K is a compact real set 
then snv zeK+m \m(z,n)\ < oo. 

Proof. We set Y = AX n for fixed n G {1, •■■ ,N}. Let z E D and be a sequence converging to z. 
Obviously exp(z p Y) — >• exp(zF) a.s. In order to conclude we need to show that the sequence (exp(z p Y j) is 
uniformly integrable. After extraction of subsequences, we can separately suppose that 

1. either min„ Re(z n ) < Re{z p ) < Re(z), for all peN, 

2. or max„ Re(z n ) > Re(z p ) > Re(z), for all p G N. 

This implies the existence of a, A G D D K such that a < Re(z p ) < A, for all p G N. 
Consequently if M > 0, for every p G N, we have 

/ — Af /> oo 

exp(yRe(z p ))d^ Y (y) + / exp(yRe(z p ))d^ Y (y) 
-oo J M 

where [iy is the distribution law of Y . Previous sum is bounded by f_ exp(ay)d/iY (y) + Jm ex P(Ay)dnY (y) 
Since M is arbitrarily big, the result is established. □ 

Lemma 3.7. Let n = 0, ■ • • ,N. 

1. E[e A - Y - - l] 2 = m(2,n) - 2m(l,n) + 1. 

2. l/ar[e AX " - 1] = m(2, n) - m(l, nf . 

3. E[e AA "- - 1] = m(l,n) - 1. 

Proof. Statements 1. and 3. follow in elementary manner using the definition of m. 

Statement 2. follows from statement 1. and the fact that E[e AXn — 1] = m(l, n) — 1. □ 

Remark 3.8. m(2, n) — n) 2 is strictly positive for any n = 1, • • • , N. In fact Assumption^ 1. implies 
that e AXn — 1 is never deterministic. 

Remark 3.9. For z G D and n G {1, • • • N}, we have E(S*) = TO ( Z > 
Proposition 3.10. for n G {1, • • • A^}, toe have 

1. AA n = E[A5„|J"„_i] = (m(l,n) - l)5„_ x . 

2. ^ar[AS n |J- n _i] = (m(2,n) - m(l,n) 2 )5 2 _ x . 
5. Condition (ND) is always satisfied. 

4- 

l — l 

" ~~ Sn-x m{2,n) - 2m(l,n) + l' 
5. The mean-variance tradeoff process K d is deterministic. 
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Proof. 1. follows from E[A5„|J"„_i] = S"„-iE[e AX ™ - 1] and Lemma O 3. 

2. Since 

E[(A5„) 2 |J- n _i] = SS-iEKe^" - l) 2 ], (3.4) 

we can write 

y ffl r[AS n |.F n _i] := E[(AS„) 2 | .F n _i] - E[A5„| J^x] 2 , 

= Sl_Me AX "-l) 2 ]-SLine AX "-l] 2 
= Sl^Var[e AX --l]. 

The conclusion follows from Lemma \'3. 71 2. 

3. We make use of Proposition 12.71 In our context we have 

E[A5 n |J-„_i] 2 (m(l,n)-l) 2 



yar[AS , n |J 7 „_i] m(2, n) — m(l, n) 2 
The denominator of the right-hand side never vanishes because of Remark 13.8 

4. It follows from (|2.3|) . (|3.4p . Lemma [3~71 1 . and point 1. of this Proposition. 

5. It is a consequence of point 3. and Definition 



(3.5) 



□ 

3.2 Discrete Follmer-Schweizer decomposition 

Similarly to |29) and |28) , we would like to obtain the discrete Follmer-Schweizer decomposition of a random 
variable of the type H = Sff, for some suitable z G C. The proposition below generalizes Lemma 2.4 of |29| . 

Proposition 3.11. Under Assumption^ let z G D fixed, such that 2Re(z) G D. Then H(z) = admits 
a discrete Follmer-Schweizer decomposition 

f H{z) n = H(z)o + E n k=iZ(z)k&S k + L(z) n 
\ H(z) N = H(z) = S z N 



where 



H{z) n = h(z,n)S z n , for all n G {0, • ■ ■ N} 

= g(z,n)h{z,n)S z l Z 1 1 , for all nG {!,••• N) (3.6) 



L{z) n = H(z) n -H(z)o-Y,£( z h&S k , for all ne{0,---N} 
and g(z,n), h(z,n) are defined by 



fc=i 



h(z,n) := ' (m(z,i)-g(z,i)[m(l,i)-l]) (3.7) 



i=n-\-l 



m(z + l,ri)-m(l,n)m(z,ri) 

9{z,n) := — — — . (3.8) 

m{2, n) — rn(l, n) z 

Remark 3.12. 1. z + 1 G D because D is convex, taking into account Assumption^ 2. 
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2. If 2Re(z) does not belong to D, for simplicity, we will set 

g(z, n) = h(z, n) = H(z) n = £(z) n = L(z)„ = 0. 

3. If K is a compact real interval, for any n G {0, • • • N} we have s\xp zGK+iR (\g(z, n)\ + \h(z, n)|) < oo. 

Remark 3.13. Suppose that (X n ) n= o,... ,jv * s a process with stationary increments i.e. such that 
X\ — Xq, ■ ■ ■ ,Xn — Xn-i are identically distributed random variables. 
According to Remark \3.4\ we have 

9M = miZ+ %~ mi «]? {Z) - (3-9) 
m(2) — rri(l) z 

We will denote in this case g[z) the right-hand side of i3.9\) . Moreover h(z,n) = h(z) N ~ n where 

h{z) = m(z) - g(z)[m{l) - 1]. (3.10) 

Proof of Provosition 1 3. TT[ Since z + 1 G D all the involved expressions are- well defined. Since L(z)o = 0, 
we need to prove the following. 

1. L(z) is a square integrable martingale. 

2. L(z)M is a martingale. 
From (|3.6p , it follows that 

AL(z) n = L(z) n - L(z) n -i = h(z,n)S* - h(z,n - 1)S^ L _ 1 - g{z,n)h(z,n)S z _ 1 (e AX ' 1 - 1); 

L(z) n is square integrable for any n G {0, • • • , TV} since 2z G D and (X n ) has independent increments. 
Since = S z n _ x e zAXn , we have 

AL(z) n = S^-i [h(z, n)e zAXn - h{z, n - 1) - g{z, n)h(z, n)(e AX " - 1)] , (3.11) 

therefore E[AL(z)„|.Fn-i] = 5^_ X E [h(z, n)e zAX " - h(z,n- 1) - g(z, n)h(z, n)(e AX " - 1)]. 

1. To show that L(z) is a martingale it is enough to show that 

E [h(z, n)e zAX " - h(z, n - 1) - g(z, n)h(z, n)(e AXn - 1)] = 0. 

Previous expression is equivalent to the relation h(z, n)m(z, n)—h(z, n—l)—g(z, n)h(z, n)(m(l, n) — 1) = 
for any < n < N which is equivalent to h(z,n — 1) = h(z,n) (m(z,ri) — g{z, n)(m(l, n) — 1)) for 
any < n < N. Previous backward relation with h(z, N) = 1 leads to (13.71) . 

2. It remains to prove that (L(z) n M n ) is a martingale. Since L(z) n and M n are square integrable for 
any n then L(z)„M n G C 1 . We prove now that K[AL(z) n AM n \J- n -i] = 0. Proposition [3TT0] 1. implies 
that the Doob decomposition 5 = M + A of S satisfies AA n = (m(l, n) — l)S n -i ■ Moreover 

AM n = AS n - AA n = S„-i(e AX " - 1) - ^(mfl.n) - 1) = £„_i(e A *» - m(l,n)). 

Coming back to (|3.11D 

AL(z) n AM n = S z n t\{e AXn -m(l,n)) 7i)e* AX ™ - n - 1) - 5 (z, n)(e AX " - 1)] . 
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Taking the conditional expectation with respect to J-'n-i, we obtain 

E[A£( z )„AM n |.F n _ 1 ] = E[5*li(e AX --m(l,n)) 

[h(z, n)e zAX " - h(z, n-l)- g(z, n)h(z, n)(e AX " - 1)] 
= S*-t\n(e AX «-m(l,n)) 

[h(z, n)e zAX " - h(z, n-l)- g(z, n)h(z, n)(e AX " - 1)] ] 
= S z n t\E[e^ AX -h(z,n) 

- e AXn h{z,n- 1) -e AX " 5 (z,n)/i(z,n)(e AX ™ - 1) 

— m(l, n)h(z, n)e zAXn + m(l, n)h(z, n — 1) 
+ ?n(l, n)g(z 1 n)h(z, n)(e AX " — 1)]. 

Again by Lemma 13.71 previous quantity equals zero if and only if 

n)m(z + 1, n) — n)h(z, n)m(2, n) — m(l, n)h(z, n)m(z, n) + m(l, n) 2 g(z, n)h(z, n) = , 

or equivalently m(z + 1, n) — g(z, n)m(2, n) — m(l, n)m(z, n) + m(l, n) 2 g(z, n) = 0. Remark 13.81 finally 
shows that g(z, n) must have the form (|3.8p . This concludes the proof of Proposition 13. Ill 



3.3 Discrete Follmer-Schweizer decomposition of special contingent claims 

We consider now options / : C —¥ K as in [55] of the type 



H = f{S N ) , with f(s) = [ s z U(dz) , (3.12) 

Jc 



where II is a (finite) complex measure in the sense of Rudin [3D, Section 6.1. An integral representation of 
some basic European calls can be found in [29] or [28] . 

The European Call option H = (St — K)+ and Put option H = (K — St)+ have a representation of the 
form p,12p provided by the lemma below. 

Lemma 3.14. Let K > 0. 

1. For arbitrary < R < 1 , s > , we have 

^ i-R+ica K l ~ z 

{S - K)+ - S= 2^l sZ zJz-^) dZ - (3 ' 13) 

^7Tl jR-ioo Z \ Z L ) 

2. For arbitrary R < 0, s > 

i pR+ioo K X ~ Z 

iK ~ s)+= 2ri sZ zJz-^) dZ - (3 ' 14) 

lm JR-ioo z \ z — L ) 

We need at this point an assumption which depends on the support of II. We set Iq := suppll n K. 
Assumption 3. 1. Ia is compact. 
2. 2I C D. 
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Remark 3.15. 1. Assumption^ is always verified (for any < R < 1 ) for the Call since Iq = {R, 1} is 
always included in [0, 1] which is a subset of by Assumption^ 2. 

2. Assumption^ is also verified for the Put, choosing suitable R provided that D contains some negative 
values. 

Remark 3.16. 1. Since D is convex, Assumption^ 2. and the fact that 2 G D imply that Iq + 1 C D. 

2. Since Iq is compact, taking II = 6 Z for some z G C, Assumption^ is equivalent to the assumptions of 
Proposition \3.11l 

3. Since Jo is compact, Assumption^ point 1. and Lemma \3.6\ imply that sup-, g2 / +iR tt.) | < oo, for 
every n = 1, • ■ • , N. 

4- Taking into account Remark \3.12\ and points 2. and 3. we also get sup zgC (|g(,z, n)| + \h{z, n)|) < oo, 
for every n = 1, • ■ • , N. 

Remark 3.17. Notice that Assumption^ is relatively weak and verified for a large class of models, whereas 
Assumption 8 required in ]28^ to derive similar results, in the continuous time setting, noticeably restricts 
the set of underlying dynamics. 

Lemma 3.18. For any n G {0, ■ • • , N}, according to the notations of Provosition \3. 1 f\ we have 

1. sup 2eC E[|ff(z)„| 2 ] < oo; 

2. sup 2eC E[|£(z)„| 2 (AS„) 2 ] < oo, for n > 1; 

3. sup 2eC E[(AL(2)„) 2 ] < oo. 

Proof. Remark 13.51 together with point 4. of Remark 13.161 show the validity of point 1. Point 3. is a 
consequence of points 1 and 2. Concerning this last point, let n G {1, • • • , N}. By Lemma [5771 1. 

E[\^ H (z) n \ 2 (AS n ) 2 } = g(z,n) 2 h(z,n) 2 E(S 2z _ 1 )(m(2,n)-2m(l,n) + l) 

= g(z, n) 2 h(z, n) 2 m(2z, n — l)(m(2, n) — 2m(l, n) + 1) 

The conclusion follows by Remark l3.16l □ 

Proposition below extends Proposition 2.5 of |29| . 

Proposition 3.19. We suppose the validity of Assumptions^ and\^ Any contingent claim H = /(S'at) 
admits the real discrete FS decomposition H given by 

f H n = H + J2 n k=1 $AS k +L% 
\ H N = H 

where 

H n = [ H(z) n U(dz) (3.15) 
Jc 

Cf = / e(*)»n(<fc) (3.16) 

Jc 

P n 

= L(z) n U(dz) = H n - H Q - YYf AS kl (3.17) 
Jc k=i 

according to the same notations as in Provosition \3.1f \ and Remark \3.1'2[ Moreover the processes (H n ) ) 
and (L^) are real-valued. 
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Proof. We proceed similarly to [29], Proposition 2.1. We need to prove that L H (resp. L H M) is a square 
integrable (resp. integrable) martingale. This will follow from Proposition 13.111 and Fubini's theorem. The 
use of Fubini's is justified by Lemma 13.181 The fact that H, £ H and L H are real processes follows from 
Remark [231 4. □ 



4 The solution of the minimization problem 

4.1 Mean- Variance Hedging 

We can now summarize the solution to the optimization problem. 

Theorem 4.1. We suppose the validity of Assumptions^ and\^ Let H = /(SW) with discrete real FS- 
decomposition 

\ H N = H. 

A solution to the optimal problem (|2.2[) is given by (Vq ,(p*) with V * = H and ip* is determined by 

<P* n = Cf + A„ ( Hn-i - H - J2 tfASi) (4.18) 



i=l 



where A„ is defined for all n 6 {1, • • • , N}, by 

\ - 1 m(l,rc)-l 
" S n _im(2,n)-2m(l J n) + l ' 1 ' 

Moreover the solution is unique (up to a null set). 

Remark 4.2. In the case that X has stationary increments, we obtain 

1 m(l) - 1 

" ~ S n - X m(2) - 2m(l) + 1 ' 

where m(n) = E(cxp(nXi)). This confirms the results of Section 2. in ]29j . 

Proof of theorem \4-l\ The existence follows from Theorem 12.111 Proposition 13.191 and Proposition 13.101 
points 3., 4. and 5. 

Uniqueness follows exactly as in the proof of Proposition 2.5 of in our case Lemma [3 . 71 gives 

Var[e AX " - 1] = m(2,n) - m(l,n) 2 . 

□ 



4.2 The Hedging Error 



The hedging error is given by Theorem 12.111 since the mean-tradeoff process is deterministic. 

Theorem 4.3. We suppose the validity of Assumptions [H and [5[ The variance of the hedging error in 
Theorem \4-l\ equals 

Jo = [ [ My,zMdy)Il(dz) , (4.20) 
Jc Jc 
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with 

j ( y z) = { s o +z ^k=iHy, z ', k )H z , k My, k )n k e=2 m (y + z ' e - *) n^U+i a C?) : y,z eauppn . ± 

: otherwise 

where 

, -s = m(2,j) -m{lj) 2 
a[J> m(2,j)-2m(l,j) + l 

and 

,/ j \ P(y,z;k)p(l,l;k) - p(y, 1; k)p(z,l; k) 

b{y,z;k) = , (4.22) 

l;fc) 

where p(y, z; fc) = to(j/ + z, fc) — k)m(z, k), y, z £ suppll. 

Remark 4.4. TTie junction p above plays an analogous role to the complex valued function with the same 
name introduced in J23j/ at Definition 4-3 in the continuous time framework. 



Proof. We proceed again similarly to the proof of theorem 2.1 of Theorem 12.111 gives that the hedging 
error is given by 

N N 

Jo = £>[(ALf ) 2 ] J] C 1 - W,') ■ (4-23) 

fc=l j=k+l 



Proposition 13.101 gives 



AA 3 = E[A5 3 -|^_i] = (m(l,j)- 1)^-1 
1 m(l,i)-l 



(4.24) 



A, = 



so 



1 - XjAAj = a(j), (4.25) 
and it remains to calculate E[(ALf ) 2 ]. bmce 



Aif = [ AL(z) k Il(dz) 
Jc 



we have 

(Aif ) 2 = / / AL(j/) fc AL(z) fe n(dy)n(dz) (4.26) 
JC Jc 

and hence by Fubini's Theorem 

E[(ALf ) 2 ] = / / E[AL(y) fc AL(z) fe ]n(^)n(dz). 

Relation (|3.1ip says that 

AL(z) fc = S»+l [h{y, k) e y AX * - h(y, k-1)- g(y, k)h(y, k)(e AX * - 1)] 
[h(z, k)e zAXk - h(z, k-1)- g(z, k)h(z, k)(e AXk - 1)] . 
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Taking the expectation we obtain 

E[AL(y) k AL(z) k ] = E[S%±f\{(h(z, k)h(y, k)m(y + z,k)- h(z, k)h(y, k - l)m(z, k) 

- h(z, k)h(y, k)g{y, k)E[e zAXk (e AXk - 1)] - h(z, k - l)h{y, k)m(y, k) 
+ h(z, k - l)h(y, k - 1) + h{z, k - l)h(y, k)g(y, k)E[e AXk - 1] 

- h(z, k)h{y, k)g(z, k)E[e yAXk (e AXk - 1)] + h(z, k)h{y, k - l)g(z, k)E[e AXk - 1] 
+ h(z, k)h(y, k)g(z, k)g(y, k)E[(e AXk - l) 2 ]}. 

Recalling that E[(e AXk - l) 2 ] = m(2, k) - 2m(l, k) + 1 and E[e AXk - 1] = m(l, k) - 1, we obtain 

E[AL(y) k AL{z) k ] = E[Sjl+*]{{h(z, k)h(y, k)m(y + z,k) - h(z, k)h(y, k - l)m{z, k) 

— h(z, k)h{y, k)g(y, k)(m(z + 1, k) — m(z, k)) — h{z, k — l)h(y, k)m(y, k) 
+ h(z, k - l)h(y, k - 1) + h(z, k - l)h(y, k)g(y, k)(m(l, k) - 1) 

(4.27) 

- h(z, k)h(y, k)g(z, n)(m(y + 1, k) - m(y, k)) 
+ h(z, k)h{y, k — l)g(z, fc)(m(l, n) — 1) 

+ h(z, k)h(y, k)g(z, k)g{y, fc)(m(2, k) - 2m(l, k) + 1)}. 

By Proposition 13. 1 il we have for x = y or z that 

h(x,k-l) = h(x,k)[m(x,k) - g(x,k)(m(l,k) -1)]. (4.28) 

We replace the right-hand sides of (|4.28|) in (|4,27[) and we factorize by h(z, k)h(y, k). Finally, after simplifi- 
cation we obtain 

E[AL(y) k AL(z) k ] = E[S^l]h{z,k)h{y,k){m(y + z,k) 

— m(z, k)m(y, k) + m(z, k)g(y, fc)m(l, k) + m(y, k)g(z, fc)m(l, k) 

— g(y, k)m(z + 1, k) — g{z, k)m(y + 1, k) 

— g(z,k)g(y, k)[m(l,k) - l] 2 

+ g(z, k)g(y, k)[m(2, k) - 2m(l, k) + 1]}. 



Hence, 



where 



and 



E[AL(y) k AL(z) k ] = E[S y k + _l\h{z, k)h(y, k)b(y, z; fc), (4.29) 

fc 

E[Sfcl*] = s v , +z E{eS y+z ' lAXk - 1 } = s^ +z Y[ m (y + z,£—l) (4.30) 



b(y, z; k) = {m(y + z,k) — m(z, k)m(y, k) — g(y, k)m(z + 1, k) — g{z, k)m(y + 1, k) 
+ m{z, k)g(y, k)m{l 1 k) + m(y, k)g(z, fc)m(l, k) 
- 9(z, k )g(y, k)m(l, fc) 2 + g(z, k)g(y, k)m(2, k)}. 

We observe that 

b(y, Z] k) = p(y, z; k) - g(y, k)p(z, 1; fc) - g(z, k)p(y, 1; fc) + g(y, k)g(z, k)p(l, 1; fc). (4.31) 
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Since, for x = y or z 

p(x,l;k) 
p(l,l;fc) 

it follows that b(y, z; k) = b(y, z; k). Finally, g^J), [[4T55) . fiTZty . ((130)) and (|OT|) give 

JV k N 

J (y,z) = sl +z Y,Kv,^k)h{z,k)h{y,k)J{m{y + z,l-l) J[ (1 - XjAAj) 

k=l 1=1 j=k+l 

N k N 

= s y Q +z Y,b(y,^k)h(z,k)h(y,k)l[m(y + z,e~l) [J a(j). 

k=l 1=1 j=k+l 

□ 

From the expression of the variance of the hedging error (|4.2ip , we can derive a sort of criterion for 
completeness for market asset pricing models. More precisely, the condition 

b(y,z;k) = 0, for all y, z € D and k G {1, • • -N} (4.32) 

characterizes the prices models that are exponential of PII for which every payoff (that can be written as an 
inverse Laplace transform) can be hedged. In the specific case of a Binomial (even inhomogeneous) model, 
we retrieve the fact that Jo(y, z) = and so Jo = 0. In fact, that model is complete. 



Proposition 4.5. Let a, b 6 R, Xj» = a with probability p^ and = b with probability (1 — pk)- Then 
Jo(y,z) = for every y,z E 

Proof. Writing p = p^, k 6 {0, 1, • ■ • , N}, we have 

p(y,z;k) = pe a{y+z) + (1 - p)e b(y+z) - (pe ay + (1 - p)e by )(pe az + (1 - p)e bz ) 
= p(l - p) ( e a{y+z) + e b{y+z) + e by+az + e bz+ay ^j 
= p(l - p) (e ay + e by ) (e az + e bz ) . 

So p(y, z; k)p(l, 1; k) = p 2 {l -p) 2 (e ay + e by ) (e az + e bz ) (e a + e b ) 2 . On the other hand, this obviously equals 
p(y,l;k)p(z,l;k). □ 

If X is a process with stationary and independent increments we reobtain the result of |29j]. 

Proposition 4.6. Let(Xk) be a process with stationary increments. We denote 



m(y) := E(cxp(yX 1 )) 

(1 

- 

m(2) - mil) 2 



m(l) — 1 

H(y) := m(y) — — ^ {m(y + 1) - m(l)m(y)) 

myz) — 171(1) 



The 



with 



m(2) - 2m(l) + 1 
Jo = [ [ My,z)U(dy)U(dz) 



J(yz) = \ ^'^ z ) B .&)=mjKr . */ a(y,z)^m( y + z) 
\ sl +z f3(y,z)Nm(y + z) N - 1 if a(y, z) = m(y + z) 



(4.33) 



16 



where 

a(y,z) = aH(y)H(z), 

m(2)m(y)m(z) — m(l)m(y + l)m(z) — m(l)m(y)m(z + 1) + m(y + l)m(z + 1) 
P\!Ji z ) — fTl\H t Z) 77-7 7Tv? ' 

m(2) — m{iy 

Proof. We observe that for k £ {0, • • • , N}, we have 

m(y + z, k) = m(y + z), h(y, k) = T-L{y) N ~ k and h(z,k) = H(z) N ~ k 

So 



n a (j) 

2=fc+l 



m(2)-m(l) J A _ a jv- fe _ 



^m(2) - 2m(l) + 1 
Consequently, expression (|4.2ip for supp(II), 



My, z) = s y +z p(y, z) £ m(y + z)^ 1 ('H{y)'H(z)a) N ~ k 
fc=i 

J+ z O/„ ~\ (m(y+z)~n(y)'H(z)a) 



Mv*)= s ° m^r-«:i if (434) 

\ s^/3(2/, z)iV™(i7 + z)^- 1 if m{y + z) = aU{y)H{z) 
This concludes the proof of the proposition. □ 



5 Numerical results 

As announced in the introduction, we will now apply the quasi-explicit formulae derived in previous sections 
to measure the impact of the choice of the rebalancing dates on the hedging error. We will consider two 
cases that motivated the present work: 

1. the underlying continuous time log-price model has stationary increments but the payoff to hedge is 
irregular, such as a Digital call, so that, as shown in [211 [24], hedging near the maturity can improve 
the hedge; 

2. the payoff is regular (e.g. classical call) but the underlying continuous time model shows a volatility 
term structure which is exponentially increasing near the maturity, such as electricity forward 
prices. For this reason it seems again judicious to hedge more frequently near the maturity, where the 
volatility accelerates. 

5.1 The case of a Digital option 

We consider the problem of hedging and pricing a Digital call, with payoff f(s) = 1\k,oo){ s ) °f maturity 
T > 0. From (35) in [22], the payoff of this option can be expressed as 

f(s) = lim — / s z dz , (5.35) 

c->oo 2m J R _ ic z 

for an arbitrary R > 0. This implies that the complex measure II is formally given by 

U(dz) = -?-^—5 R (dRe(z))d(i(Im(z))). (5.36) 

2lTl z 

However, such measure is only cr-finite so that application of Theorem l4.1l is not rigorously valid. Nevertheless, 
using improper integrals one is able to recover an exploitable form for applications. 
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Proposition 5.1. Let f(s) = 1[k,oo[{s)- We suppose again the validity of Assumption^ and\3i Let R > 0. 

1. The FS-decomposition of the contingent claim H = /(Sn) is given by 

\ H N = H 

where 

H n = lim / H{z) n n{dz) (5.37) 



(— J-oo 



R+i[-i 



£—>oo 



£ - lim / f(«)„n(<fo) (5.38) 



~ n 

L* = L(z) n U(dz) = H n -H -y2^AS k , (5.39) 

Jc k=i 



according to the same notations as in Provosition \3.1f\ and Remark \3.12\ 

2. The solution to the minimization problem is still given by Theorem \4-l\ 

3. The variance of the hedging error is given by 

lim / / J (y,z)Ile(dy)ILe(dz), 
Jc Jc 

where for each £ > ; Hi is the finite complex measure defined by YLi{B) = H(B f] (iJ + i[—£, £])) for a 
Borel set BcC. 

Proof. We proceed similarly as in Lemma 4.2 of [55] . For £ > 0, we denote fn : C — > C denned by 
fi(s) := J c s z Hi(dz). According to the proof of Lemma 4.2 of (55], there is it G R such that 

\f{s)-h{s)\ <us R , VseR 

The Lebesgue's dominated convergence theorem implies that lim^oo E [(/(.(Sn) — /(£W)) 2 ] = 0. Setting 
H e = fe(S N ), H = f{S N ) we get limbec E [(H - H 1 ) 2 ] = 0. Item 1. follows by Proposition O 
Item 2. follows by the same arguments as the proof of Theorem 14. II 

Item 3. follows exactly as in step 3. of the proof of Lemma 4.2 of [215] . □ 

In this section, this will be assumed so that formula (|4.20|) will be used in the case of a Digital option. 
The underlying process S c is given as the exponential of a Normal Inverse Gaussian Levy process (see 
Appendix H B. i.e. for all t £ [0,T], 

S% = e Xt , where X c is a Levy process with XI ~ NIG(a, (3, <5, fi) . 

Given N + 1 discrete dates = to < t\ <•••<%= T, we associate the discrete model pricing X = X N 
where Xk = A t c fc , k £ {0, . . . , N}. A is a discrete time process with independent increments. The related 
cumulant generating function z i-> m{z,k) associated to the increment AX)- = X k — X k —i = X£ k — X£ 
for k £ {1, • • ■ N} is defined on D = [—a — (3;a — /?]. We refer for this to [55] Remark 3.21 2., since X c is a 
NIG process. By additivity we can show that 

m(z, k) = E[cxp(zAX fe )] = cxp (kt k [/j,z + S( a 2 - (3 2 - ^a 2 - (f3 + z) 2 )} \ , for z £ D, k £ {0, . . . , N}. 

(5.40) 
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For other informations on the NIG law, the reader can refer to Appendix [6] B. 

Assumption [5] 1. is trivially verified, Assumption [5] 2. is verified as soon as 2 < a—0. Thanks to Remark l3.15l 
Assumption [3] is automatically verified for the Call and Put representations given by Lemma 13.141 and, by 
similar arguments, even for the digital option. 

The time unit is the year and the interest rate is zero in all our tests. The initial value of the underlying 
is so = 100 Euros. The maturity of the option is T = 0.25 i.e. three months from now. Four different 
sets of parameters for the NIG distribution have been considered, going from the case of almost Gaussian 
returns corresponding to standard equities, to the case of highly non Gaussian returns. The standard set of 
parameters is estimated on the Month-ahead base forward prices of the French Power market in 2007: 

a = 38.46 , = -3.85 , 6 = 6.40 , fj, = 0.64 . (5.41) 

Those parameters imply a zero mean, a standard deviation of 41%, a skewness (measuring the asymmetry) 
of —0.02 and an excess kurtosis (measuring the fatness of the tails) of 0.01. The other sets of parameters 
are obtained by multiplying parameter a by a coefficient C, (J3, 8, fi) being such that the first three moments 
are unchanged. Note that when C grows to infinity the tails of the NIG distribution get closer to the tails of 
the Gaussian distribution. For instance, Table [T] shows how the excess kurtosis (which is zero for a Gaussian 
distribution) is modified with the four values of C chosen in our tests. We compute the Variance Optimal 



Coefficient 


C = 0.14 


C = 0.2 


C = \ 


C = 2 


a 


5.38 


7.69 


38.46 


76.92 


Excess kurtosis 


0.61 


0.30 


0.01 


4. 10~ 3 



Table 1: Excess kurtosis of X\ for different values of a, (j3, 5, fj.) insuring the same three first moments. 

(VO) hedging error given by (|4.20[) . for different grids of rebalancing dates. The corresponding initial capital 
Vq denoted by V * = Ho in Theorem 14. II is computed using Proposition ^. 191 

In particular, we consider the parametric grid introduced in [21] and [24] ir b ' N := {0 = t b ' N , t b { N , • • • , t b ^ N } 
defining, for any real b E (0, 1], N rebalancing dates such that 

t b,N = T _ t q _ k y /b for all fc e { 0j ... ;7V _ !} , ( 5 .42) 

Note that tt 1 ^ coincides with equidistant rebalancing dates whereas when b converges to zero, the rebalancing 
dates concentrate near the maturity. To visualize the impact of parameter b on the rebalancing dates grid, 
we have reported on Figure Q] the sequences of rebalancing dates generated by -K h N for different values of 6. 
We have reported on Table [2] the standard deviation of the Variance Optimal hedging error for different 
values of coefficient C and different choices of rebalancing grids. More precisely, we have considered three 
types of rebalancing grids, for N = 12 rebalancing dates. 

1. Equidistant rebalancing dates (corresponding to ir 1,N ); 

2. TT b ' where b* is obtained by minimizing the Variance Optimal hedging error w.r.t. to parameter b; 

3. The non parametric optimal grid n* obtained by minimizing the Variance Optimal hedging error w.r.t. 
the N rebalancing dates. 

Notice that in both cases the optimal (parametric and non parametric) grid is estimated by an optimization 
algorithm based on Newton's method. 
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Figure 1: Sequences of rebalancing dates for different values of 6, for N = 12. 



First, one can notice that for any choice of rebalancing grid, the hedging error increases when C decreases. 
Hence, one can conclude, as expected, that the degree of incompleteness increases when the tails of log- 
returns distribution get heavier. 

Besides, one can notice that the parametrization (|5.42[) of the rebalancing grid seems remarkably relevant 
since the optimal parametric grid 7r & achieves similar performances as the optimal non-parametric grid it* . 
Moreover, we observe that the hedging error can be noticeably reduced by optimizing the rebalancing dates 
essentially for C > 1 i.e. around the Gaussian case. In these cases, one can observe on Figure [5] that 
the optimal rebalancing grid is noticeably different from the uniform grid since rebalancing dates are much 
more concentrated near maturity. This confirms the result of |21| that shows that, in the Gaussian case, 
taking a non uniform rebalancing grid (corresponding to b — 0.5) allows to obtain a hedging error with 
the convergence order for the L 2 norm of TV" 1 / 2 (up to a log factor) improving the rate iV -1 / 4 achieved 
with a uniform rebalancing grid (i.e. 6 = 1), obtained in [26] . However, it is interesting to notice that this 
phenomenon is less pronounced when the tails of the log-returns distribution get heavier. In particular, one 
can observe on Figure [3] that the hedging error gets less sensitive to the rebalancing grid when C decreases 
even if the optimal grid seems to get closer to the uniform grid. 

5.2 The case of electricity forward prices 

We consider the problem of hedging and pricing a European call, with payoff (F^ d — K)+, on an electricity 
forward, with a maturity T = 0.25 of three month. The maturity T is supposed to be equal to the delivery 
date of the forward contract T = Td- Because of non-storability of electricity, the hedging instrument is the 
corresponding forward contract. Then we set Sf = F^F , where the forward price F T is supposed to follow 
the NIG one factor model (jl.ip with m = 0, ox = and a s — a > 0. This gives 



S£ = e x °, where X£ = / ae~ x{T ~ u) dA u where A is a NIG process with A x ~ NIG(a, j8, S, fi) . 

Jo 

(5.43) 
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O = I 


u = 1 


/~» no 


O = U.14 


10 x STDvq^*) 


1.483 (30.82) 


1.652 (34.33) 


2.663 (54.80) 


3.017 (61.53) 


10 x STDvo^,) 


1.520 (31.58) 


1.685 (35.01) 


2.665 (54.84) 


3.017 (61.53) 


10 x STDvow 


1.892 (39.32) 


1.952 (40.56) 


2.691 (55.38) 


3.028 (61.76) 




0.4903 


0.4859 


0.4813 


0.4812 


V (tt*) 


0.4903 


0.4860 


0.4814 


0.4813 


b* 


0.4078 


0.4394 


0.6106 


0.6710 



Table 2: Standard deviation of the Variance Optimal hedging error (xlO) (reported within parenthesis in percent 
of the initial capital Vb^ 1 )), initial capitals for b —\ and 6 = b* , optimal grid parameters b* , for different choices of 
parameters C with N — 12 and K = 99 (Digital option). 




Figure 2: Parametric and non parametric optimal rebalancing grids for different choices of parameter C with iV = 12 
and K = 99 (Digital option). 
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Figure 3: Standard deviation of the Variance Optimal hedging error as a function of b, for different choices of 
parameter C (b* being indicated by the dashed line abscissa) with N = 12 and K = 99 (Digital option). 

Given N + 1 discrete dates = to < t\ < ■ ■ ■ < tjy = T, we consider the discrete process X = X N where 
Xk = Xt k ,0 < k < N. We denote again by z H m(z,k) the cumulant generating function associated with 
the increment AXk = Xk — Xk-i for k G {!,•• -N}. That function and its domain can be deduced from 



Lemma 3.24 and Proposition 6.2 in [25], see also ([6"37|) . The domain D contains D := [-^, ^] 
and given for any z € D, k = 0, . . . , N, 



cxp 




m(z,k) = E[exp(z / ae~ x{T ~ u} dA u )} 

c-l 

K A (z u du I , with z u = zcre" A(T_u) 



exp / [/jlZu + S( y/a 2 - P 2 - V® 2 ~W + z u ) 2 )]du , (5.44) 



where k a is recalled in formula (|6.57p . Hence Assumption [2J 1. is obviously satisfied since A ^ and 
Assumption [2J 2. is verified as soon as a < thanks to Remark 13.151 Assumption [5] is automatically 

verified for the call representation given by Lemma 13.141 

Parameters are estimated on the same data as in the previous section, with Month-ahead base forward prices 
of the French Power market in 2007. For the distribution of Ai this yields the following parameters 

a = 15.81 , P = -1.581 , 5 = 15.57 , fi = 1.56 , 

corresponding to a standard and centered NIG distribution with a skewness of —0.019 and excess kurtosis 
0.013. The estimated annual short-term volatility and mean-reverting rate are a = 57.47% and A = 3. 
We have reported on Figure the standard deviation of the hedging error as a function of the number of 
rebalancing dates for four types of hedging strategies. 

• Variance Optimal strategy (VO) with the uniform rebalancing grid (dark line) and with the op- 
timal rebalancing grid ir* (dark dashed line). Both variances are computed using formula (j4.20[) 
applied to the process (|5.43[k 
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• Black-Scholes strategy (BS) implemented at the discrete instants of the uniform rebalancing grid 
(light line) and of the rebalancing grid tt* (optimal for the Variance Optimal strategy) (light 
dashed line). Both variances are computed using Theorem 3.1 of [3J extended to non-stationary log- 
returns, to derive a quasi- explicit formula for the variance of the BS hedging error. Indeed, in [2J, the 
authors uses the Laplace transform approach, to derive quasi-explicit formulae for the mean squared 
hedging error of various discrete time hedging strategies including Black-Scholes delta when applied to 
Levy log-returns models. This extension of this result to the general case when X is a non-stationary 
process with independent increments is given below. 

Proposition 5.2. Let v be an admissible strategy satisfying 

r(z) n S z n z\Tl(dz) (5.45) 

for n = 1,...,N, where f v {z) n is a deterministic function of the complex variable z. Let c be the 

ini, 

by 



initial capital; the bias and the variance of the hedging error e(v, c) := H — c — Ylk=i v kASk is given 



E[e(t/,c)] = / S * 



N N 



J m(z, k)-J2 f V (z)k(m(l, k) - 1) JJ "*(*, 1 - K 



IL{dz) - c (5.46) 



_fc=l fe=l 1=2 

E[e(u,0) 2 ] = J J S y Q +z (v 1 (y,z)-v 2 (y,z)-v 3 (y 1 z) + v 4 (y 1 z))n(dz)Il(d y ) , (5.47) 

where 

N 

vi(y, z ) = Y[m(y + z,k) 

k=l 

N k-1 N 

V2(y,z) = ^ f v {y)k [m(z + 1, fc) - m(z, k)) JJ m(y + z, I) JJ m(z,l) 

fe=l (=1 l=k+l 

N k-1 N 

va(y,z) = ^2f v (z)k[m(y + l,k)-rn(y,k)][[rn(y + z,l) JJ m(y,l) 

k=l 1=1 l=k+l 

N k-1 

v 4 (y,z) = J2f V (y)kf v (z)k[m(2,k)~2m(l,k) + l}l[m(y + zJ) 

k=l 1=1 

N k-1 j-1 

+ E E f v {z)kf v {y) 3 [m(y + 1, fc) - m(y, fc)] JJ m(y + z,l) JJ m(y, l)(m(l,j) - 1) 

j=2 k<j 1 = 1 l=k+l 

N j-1 k-1 

+ EE^W^^ W« + 1.J) - ^.i)] II m (tf + II rn(z,l)(m(l,k)-l) . 

3=2 k>j 1 = 1 1=3+1 

Therefore, the variance of the hedging error is 

Var(e(w; c)) = Var(e(w; 0)) = E[(e(u; 0) 2 ] + E[e(v; 0)] 2 . 

Proof. The proof is similar to the one of Theorem 3.1 of [2]. □ 
Remark 5.3. In the case of Black-Scholes delta hedging strategy 

N 



A*)« = - IJ k ) ' where "»*"(*. k ) = CX P (- Var f Xk] z + Var ^ Xk K 

k=n ^ 
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Observing Figure [4] one can notice that, as expected, in all cases, the hedging error decreases when the 
number of trading dates increases. Observing the continuous lines, corresponding to a uniform rebalancing 
grid, one can notice the remarkable robustness of the Black-Scholes strategy. Indeed, in spite of the non 
Gaussianity of log-returns and the discreteness of the rebalancing grid, the Black-Scholes strategy is still 
quasi optimal in terms of variance. 

Besides, in this case, the impact of the choice of the rebalancing grid seems to be more important than 
the choice of log- returns distribution (Gaussian or Normal Inverse Gaussian). For instance, using the VO 
strategy with the optimal rebalancing grid ir* instead of ir 1 allows to reduce 9% (for TV = 10) of the hedging 
error standard deviation. The BS strategy shows similar performances to the VO case, when implemented 
at the rebalancing times tt* . Indeed BS optimal rebalancing grid (in terms of variance) appears to be close 
to 7r* (up to 10~ 4 ). Moreover, one can observe on Table [3] that here again, the parametrization (|5.42|) of the 
rebalancing grid seems to be particularly well suited since it achieves minimal hedging errors comparable to 
the one achieved with the nonparametric optimal grid n* . 

Notice that our analysis only considers the variance of the hedging error. To obtain the mean square 
error, one should add the bias contribution which is of course zero for the variance optimal strategy but it 
is in general non negligible for the Black-Scholes strategy. In particular, we can observe that this bias term 
varies strongly with the parameters of the NIG distribution. 

For instance, for TV = 2 uniform rebalancing dates, replacing parameter /3 by —/3 increases the bias 
(defined as (|5.46[> . with initial capital c = V® 8 ) from -0.04 to 4.45. Moreover, one should also observe 
that the drift and the skewness of log-returns also impact the standard deviation of the BS hedging error. 
Changing again /? by — j3 implies an increase of the log-returns expectation (resp. skewness) from to 3.12 
(resp. from -0.02 to 0.02) which induces an increase of the standard deviation of the BS hedging error from 
4.91 to 5.92, whereas the standard deviation of the VO hedging error decreases from 4.83 to 2.10. 




1 - 

0.5- 

nl i i I i i i i i i I 

5 10 15 20 25 30 35 40 45 50 

Number of trading dates 

Figure 4: Standard deviation of the hedging errors as a function of the number of rebalancing dates TV, for K = 99 
(Call option). 

To analyze the impact of the rate of volatility increase on the optimal rebalancing grid, we have computed 
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N = 2 


N = 5 


N = 10 


N = 25 


N = 50 


STDyQ(^*\ 


4.5683 (53.23) 


3.1129 (36.10) 


2.3807 (27.56) 


1.7790 (20.57) 


1.5233 (17.61) 


STDyo^b*) 


4.57167 (53.27) 


3.1550 (36.59) 


2.4186 (28.00) 


1.8023 (20.84) 


1.5354 (17.75) 


STD VO (iri) 


4.8331 (56.32) 


3.4012 (39.44) 


2.6154 (30.28) 


1.9275 (22.29) 


1.6145 (18.66) 




4.9137 (57.26) 


3.4196 (39.66) 


2.6217 (30.35) 


1.9329 (22.35) 


1.6231 (18.76) 




4.6291 (53.94) 


3.1273 (36.27) 


2.3884 (27.65) 


1.7886 (20.68) 


1.5344 (17.74) 




8.5818 


8.6232 


8.6380 


8.6469 


8.6499 


V (n*) 


8.5895 


8.6275 


8.6406 


8.6493 


8.6531 


b* 


0.5917 


0.6298 


0.6284 


0.6203 


0.6172 



Table 3: Standard deviation of the Variance Optimal hedging error (reported within parenthesis in percent of the 
initial capital Vo(ir )), initial capitals, optimal grid parameters b* , for different choices of rebalancing dates N (Call 
option). 

the hedging error standard deviation for several values of parameter A choosing the corresponding volatility 
parameter a such that Var^Xx) = fr(l — e~ 2AT ) is fixed. The resulting pairs (A, a) are reported on Table[4] 
Coupling those parameters allows us to obtain comparable options for different parameters A; at least this 
ensures a fixed initial capital in the BS framework (with Vq = 8.7037). On Figure we have reported the 
optimal grid parameter b* minimizing the standard deviation of the VO hedging error for different values of 
A. As expected, when A increases, i.e. when the volatility increases more rapidly near the maturity, then b* 
decreases indicating that the optimal rebalancing dates concentrate near the maturity. On Figure one can 
observe that the hedging error increases with A even when the rebalancing dates are optimized. However, 
optimizing the rebalancing dates allows to reduce noticeably the hedging error, specifically for high values of 
A. For instance, it allows to reduce 7.5% of the error standard deviation when A = 3 and 17.9% when A = 9. 



A 


1 


2 


3 


6 


9 


a 


0.4662 


0.5202 


0.5747 


0.7349 


0.8823 



Table 4: Short term volatility a (s.t. Var(XT) = fr(l — e 2XT ) is fixed) for different values of parameter A with 
N = 10 and K = 99 (Call option). 
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A. (year 1 ) 



Figure 5: Optimal rebalancing grid parameter b* as a function of A for K = 99 and N = 10 (Call option). 




Figure 6: Standard deviation of the hedging error as a function of A for K — 99 and N = 10 (Call option). 
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6 Appendix 



A: A general convergence theorem for FS decompositions 
Proposition 6.1. Let (-ff )^gNu{+oo} ^ e a se Q uen ce of r.v. in C 2 (fl, J~n)- Let 

1 K = H* 
be the FS- decomposition of H e . Suppose that H e — > H°° in C 2 (il). Then, for £ — > +oo. 

1. Hq — > mC 2 (T Q ); 

2. — > in probability for any n G {1, . . . , iV}; 

3. L e N -> L™ in C 2 (n). 

Proof. For nE {!,..., N}, £ eNU {+00} we have 



(6.48) 



Hi = H 



CAS n + ALi 



(6.49) 



For technical reasons we set := and L £ N+1 := L l N . The result will follow if for every n E {0, . . . , N}, 

for £ — > +00 we have 



1. m -> H™ in £ 2 , 



2. E 



We will prove l.,2. and 3. by backward induction on n e {0, . . . , TV} starting from n = N. The step TV 
of the induction is constituted by the assumption, in particular 1. and 3. are verified by assumption and 2. 
is trivially verified. 

Suppose that l.,2. and 3. hold for some n E {1, . . . , N}, we will prove their validity for the integer n — 1. 
First, 1. implies that E [i^lJ^i] — >e_> +00 E [H^\F n ^i\ in £ 2 (Q). We continue taking the conditional 
expectation with respect to F n -i m (|6.49p . This gives 



^[Hi\F n - l }=H l n _ 1 +CAA n . 
The difference between (|6.49[) and (|6.50p gives 

H* - E \Hi\Tn-!] = CAM n + ALi, * e N U {+^}. 

Consequently 



(6.50) 



So 



E 



< -H™=E [Hi ff-lJn-r] = - ft ) Atf„ + A(Ll - 
« - ff„°°) -E [i# - ff~|^„-i] 2 l = E [(£ - ^T 3 ) 2 (AM n ) 2 l + E [A (L e n - L~)] 2 : 



(6.51) 



in fact 



E ((£ - C) AM„A (l£ - L~)) = E ((£ - £°)E ((AM„)A(i£ - L~)|^ n _i)) = 0, 
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because M. (l/ — is a martingale. Since the left-hand side of ()6.51|) converges to zero when I — > oo, it 
follows that 



E 



{AM n ) 2 
E[A «-L r T)] 



(6.52) 



This shows 2. and 3. of the (n — l)-step of the backward induction. It remains to show item 1. By (|6.49D . 
we have 

H n-l ~ H n-1 = H n - H ™ - &S n (£ n - - A (L„ - L£°) . 

Since _ff r £ t — and A (L^ — L^°) converge to zero in C 2 , it remains to show that AS n (£„ — £^°) — >£^oo 
in C 2 (Vl) when £ — > oo. Now AM n — — S^-j-oo in C 2 (Vt) and so by (|6.52[) we only have to prove that 



E 



0. 



By the (ND) condition and item 1. of Remark 12.41 we have 

(A A n ) 2 = (E(A5„|J-„_ 1 )) 2 <<5E((A5„) 2 |J-„_ 1 



Consequently 



= S[(AA n ) 



(AA n ) 2 < 



1-5 



3 



(AM„) |J"„_i 



(AM„) |J"„_ 



So the left-hand side of (|6.53[) is bounded by 

6 



1 - S 



:3 



(AM„)' 



The result is finally established. 



(6.53) 



□ 



B: The Normal Inverse Gaussian distribution 

The Normal Inverse Gaussian (NIG) distribution is a specific subclass of the Generalized Hyperbolic family 
introduced by Barndorff-Nielsen in 1977, see for instance [3]. The density of a Normal Inverse Gaussian 
distribution of parameters (a, j3, S, /i) is given by 



f NIG (x-,a,(3,S,^ = -exp(6^a 2 -(3 2 +(3(x-fi)) Kl ^ ad ^l^ {x =jff^ , for any x e R , (6.54) 

7T V 1 + (' X - M) /<•> 



where i\i denotes the Bessel function of the third type with index 1 and where the parameters are such that 
5 > 0, a > and a > |/3|. Afterwards, NIG(a, (3, S, n) will denote the Normal Inverse Gaussian distribution 
of parameters (a, /?, S, fx). 

A useful property of the NIG distribution is its stability under convolution i.e. 

NIG(a, f3,5i, fti)* NIG(a, (3, 8 2 ,/i 2 ) = NIG(a, p, <5i + 5 2 , m + fx 2 ) . 

This property shared with the Gaussian distribution allows to simplifies many computations. 

If X is a NIG (a, f3,S, fi) random variable then for any a £ K + and b € K , Y = aX + b is also a NIG 
random variable with parameters (a/a, /3/a, a5, afx + b). 
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The mean and the variance associated to a NIG(a, ft, 5, fx) random variable X are given by, 

EX = /i + — , VarV = , with 7 = ^/ a 2 - /3 2 . (6.55) 
The characteristic function of the NIG distribution is given by cxp^ nig) where ^nig verifies 

*jvre(ii) = logE[exp(raJT)] = ifxu + 5{s/a 2 - P 2 - sJol 2 - (/3 + iu) 2 ) for any ueR. (6.56) 
The moment generating function of the NIG distribution is particularly simple, 

k a (z) = k nig (z) = logE[exp(zA)] = fxz + S^a 2 - P 2 - ^/a 2 - (/? + z) 2 ) , for Re(z) 6 [-(« + /?); a - /?] . 

(6.57) 

The Levy measure of the NIG distribution is given by 

Sex. 

F NIG (dx) = e l3x —-K 1 {a\x\) dx for any iel. (6.58) 
7r|x| 

Notice that the Levy measure does not depend on parameter fx. 
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